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Abstract 

We propose a new parametrization of the B — > ir vector form factor, f + (q 2 ), as an expansion in 
powers of a conformal mapping variable, which satisfies unitarity, analyticity and perturbative QCD 
scaling. The unitarity constraint is used also for defining the systematic error of the expansion. We 
fit with the new parametrization the available experimental and theoretical information on exclusive 
B — > nlis decays, making a conservative estimate of the effects of correlations in the systematic and 
statistical errors of the lattice results. With four parameters to describe fj r {q 2 ), the systematic 
error is negligible in the whole semileptonic region. We also obtain \V u b\ = (3.50 ± 0.24) x 10 -3 
where, in our approach, the uncertainty is predominantly statistical. 

PACS numbers: 12.15.Hh, 13.20.He 
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I. INTRODUCTION 



As shown recently 



the extraction of \V u b\ from the exclusive semileptonic 
B — ► 7r decays has become competitive with determination from inclusive decays. The 
exclusive decay approach requires the theoretical description of the matrix element 



{'K{p^)\u^ IJ b\B{p^ + q)) = 2p n ^ + q^ - q^ 



m 



/ + (<? 2 ) + v 



/o(g 2 



i) 



where g is the momentum of the lepton pair and f+(q 2 ) and foil 2 ) denote the vector and 
scalar form factors, respectively. For light leptons, I = e, fi, only the vector form factor 
contributes to the spectrum 
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A 3/ V)l/ + (s 2 )P 



(2) 



B 



where A(g 2 ) = (m 2 B + m 2 — q 2 ) 2 — 4m 2 3 m 2 is 4m 2 B times the pion three-momentum squared 
in the 5-meson rest frame. The physical range of semileptonic decays is < q 2 < t_, with 
t_ = (m 2 B0 - m w+ ) 2 = 26.42 GeV 2 . 



The Q 2 -spectrum of B — > -niv decays has been measured with increasing precision by 

inn nn 

the CLEO [5|, |6j, Belle [7[, and BaBar collaborations [8|, |9|. The accuracy of theoretical 
calculations of the form factors is also continuously improving. The calculations are based 
either on QCD light-cone sum-rules (LCSR), which provide reliable determinations at small 
12| . or on lattice simulations, which give accurate results at larger values of q 2 
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From general principles of quantum field theory it is known that the vector form factor 
f+(q 2 ) is a real analytic function in the complex g 2 -plane cut for q 2 > t + , where y/t^ = 
(rriBO + m n +) = 5.419 GeV is the Bit threshold. Angular momentum conservation imposes 
the behavior Im/ + (g 2 ) ~ (q 2 — t + ) 3 ^ 2 near the threshold. In addition, f+(q 2 ) has a pole 
below the branch point, at q 2 = m 2 Bt {m B * = 5.325 GeV). As shown in [19|, unitarity applied 
to a certain QCD correlator provides a constraint on the magnitude of f+(q 2 ) along the cut, 
while scaling in perturbative QCD requires a power-law fall-off, like 1/q 2 , up to logarithmic 



corrections 
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2l|. 



To maximize the usefulness of this information about B — > ttEu decays and, for instance, 
determine \V u b\, it is crucial to have a parametrization of f+{q 2 ), which satisfies the above 
theoretical requirements, and whose associated systematic error is quantifiable and small 
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compared to experimental and other theoretical errors. The purpose of the present work 
is to provide such a parametrization. In Sec. [Til we discuss the recent models proposed in 
the literature, and show that they do not satisfy all of the constraints mentioned above. In 
Sec. IIII[ we propose a simple analytic parametrization of f + (q 2 ), which combines the pole 
factorization with an expansion in powers of a conformal mapping variable, and in Sec. IIVI 
we express the unitarity constraints in terms of the coefficients of this expansion. By fitting 
the experimental differential decay rate B — > irlv and the values of f+{q 2 ) calculated from 
LCSR and the lattice, we obtain in Sec. IVIIII a model-independent representation of the 
form factor f+(q 2 ) in the physical region. The procedure also yields a precise value of \V u b\, 
given in Sec. [DO 



II. RECENT PARAMETRIZATIONS OF THE FORM FACTOR 



A comparison of the various parametrizations proposed in the literature was done recently 
in |22l| . As shown there, the simple expressions with a pole at q 2 = m 2 B *, proposed in 
[ill . Q], cannot provide an accurate description of the form factor in the whole physical 
domain. Thus, more systematic expansions, which incorporate the constraints of analyticity 
and unitarity, have been considered. 

A first type of parametrization is inspired by the technique of unitarity bounds proposed 



by Okubo [2J, |25j and applied subsequently to various semileptonic form factors 
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281 ]. They are obtained by exploiting the analyticity and positivity properties of vacuum 
polarization functions of the type U^ u (q)=i J d 4 xe iq ' x (0\T{Jjx)Jl(0)}\0) , where J M = u^b 



here. In this approach, the form factor can be expanded as 

1 



29] 



U(q 2 



5(g 2 )0(g 2 ,t o ) 4„ 
where a n (t ) are real coefficients and z = z(q 2 ,t ) is the function 

z(q 2 ,t ) 



(3) 



Vu - <? 2 - 



<t + - g 2 + ^/U=T Q 



(4) 



which maps the g 2 -plane cut for q 2 > t + onto the disk \z(q 2 , to)] < 1 in the z complex plane, 
such that z(t + ,t ) = —1 and 2(00, t ) = 1. The parameter t < t+, which is arbitrary, 
determines the point q 2 mapped onto the origin in the z plane, i.e. z(to, to) = 0. 



3 



The function B(q 2 ) is the Blaschke factor 

B(a 2 ) - *(gVo) -*( m g-»*o) _ ( 2 2 \ /,n 

which accounts for the pole at g 2 = m 2 B , . By construction |5(g 2 )| = 1 for q 2 > t + . 
The last factor in (j3J), the outer function (f>(q 2 ,t ), has the expression 



1 t — a 2 



32tt Xi -(0) vv + * v + u '(t+-<o)i/< 



x( V / ^ r 7+ v 7 ^)- 5 (V / V r 7+ v 7 ^ 3 ^) 372 , (6) 



where Xi-(O) i s the derivative of the transverse component of the polarization function 



n M i,(<?) at the Euclidean momentum Q 2 = —q 2 = [191 ] . Because of the large value of 



the 6-quark mass, this quantity can be computed by means of perturbative QCD and the 



operator product expansion [19j, [30|. On the other hand, the spectral function associated 
with Il liv {q) is a sum of positive contributions. Thus, if we assume that it is saturated by Bit 
intermediate states, unitarity and crossing symmetry guarantee that the coefficients a n (to) 
satisfy the inequality 

CO 

E°n(*0)<l. (7) 

n=0 

This allows one to calculate bounds on the values of the form factor or its derivatives at 



points inside the analyticity domain, in particular in the physical region. For the B — > tt 
form factors, such bounds were investigated in [l| and jl9j ]. 



The expression ()3]) was also adopted as a parametrization of the form factor i n |22l . l3ll . [32] . 
In this case the expansion is truncated at a finite order. However, as noticed in {33], the form 
factor increases then like f+(q 2 ) ~ (q 2 ) 1 ^ 4 at large \q 2 \, in contradiction with perturbative 
QCD scaling. This behavior follows from the expression ([6]) for the outer function, taking 
into account that all the other factors in Eq. ([3]) are finite at \q 2 \ ^00. Moreover, when the 
series is truncated, the expression has an unphysical singularity at the Bti production 
threshold t + , produced by the factor (t + — q 2 ) in the numerator of Eq. ([6]). This unphysical 
singularity may distort the behavior near the upper end of the physical region, where the 
form factor is poorly known. 

It should be noted that in the calculation of bounds one uses the full expansion in Eq. (JSj), 
with an infinite number of terms. Then the series cancels the zeros of the function 0(g 2 ,to) 
at q 2 = t + and at q 2 — > 00, restoring the required properties of f + (q 2 ). Actually, it can be 



shown that imposing the condition that the series vanishes at threshold or at infinity does 
not change the unitarity bounds (this answers a question raised in [33| about the possibility 
of improving the bounds in this way). 

A second type of parametrization, used recently for the B — > tx form factors in [2, 4 1, is 
based on the Omnes representation 34J , which expresses an analytic function in terms of its 
phase along the boundary of the analyticity domain. If the phase 5(t), defined by f+(t+ie) = 
| exp(id(t)) for t > t + , has a finite limit at infinity, the Omnes representation requires 
only one subtraction. Taking into account the pole at q 2 = m 2 B * and assuming, as in 
that the form factor does not have zeros in the complex plane, the representation reads 



a, 



(? 2 



2 2 
/ 2\ m B* ~ Ql 
-(Qi)— 2" eX P 

m B* ~ 1 



q 2 -q\ f 8(t) dt 
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(t-q 2 )(t-q 2 



(8) 



where q\ is an arbitrary subtraction point. By Watson's theorem [35[, the phase 5(t) is 
equal, below the first inelastic threshold, to the phase of the P-wave with I = 1/2 of the 
ttB — > ttB elastic scattering. Since this phase is not known, in Refs. [2J, |4J] the contribution 
of the integral was suppressed by using a multiply-subtracted dispersion relation. Neglecting 
altogether the dispersion integral, the form factor is represented in ^] as 

1 n 

iW) = — — - 2 H[f+(®( m B* -<zf)r (<?2) 

171 B* ~ 

where 



Q 1 



(9) 



3=1 



«i(9 2 )= II 



r - Qt 



Qt 



(10) 



However, it is easy to see that the expression ([9]) defines an entire function in the complex 
g 2 -plane, with no cut for q 2 > t + . So, this parametrization does not have the proper 
structure of the physical form factor required by analyticity and unitarity. Moreover, the 
expression (Q exhibits at \q 2 \ —>■ oo an exponential behavior like exp[C(g 2 ) n ~ 1 ], where C 
is a combination of the values f + (q 2 ). This anomalous behavior follows from the multiple 
subtraction of a dispersion relation that requires in general only one subtraction. The values 
f+(q 2 ) for j > 1 are not independent: according to (jBJ), they can all be expressed in terms of 
f+{qf) and the dispersion integral. By taking into account these relations in the multiply- 
subtracted dispersion relation, one recovers the original relation (jBJ). However, if the values 
f+{q 2 ) are treated as independent, the form factor behaves as an exponential at large q 2 , in 
contradiction with QCD scaling. 



Though the shortcomings discussed in this section formally concern the behavior of the 
parametrizations of Eqs. ([3]) and (Q outside the semileptonic domain, it is important to 
construct a representation of the form factor which has the correct analyticity properties in 
the whole complex plane. Besides the obvious argument that a parametrization that does 
not satisfy these properties cannot describe the form factor correctly, the introduction of 
unpysical singularities, sometimes close to the semileptonic domain, can distort the form 
factor in that region. Given the levels of precision currently reached and expected in the 
study of exclusive B — > n decays, such distortions are unacceptable. 



III. A NEW PARAMETRIZATION FOR f + (q 2 ) 

We start by the remark that the product (1 — q 2 /m 2 B *)f + (q 2 ) is analytic in the complex 
g 2 -plane cut along the real axis for q 2 > t + and is finite for q 2 —>■ oo, due to the scaling 
behavior f+(q 2 ) ~ l/? 2 - An expansion of the product that converges in the whole complex 

Qne is obtained in terms of a variable that conformally maps the cut q 2 plane onto a disk 
The variable z = z(q 2 ,to) defined in Eq. (TJJ performs precisely this mapping. Thus, 
we propose the simple parametrization 

utf) = l i /m2 X>(*o)* fc . (ii) 

The polynomial in powers of z displays the branch point at q 2 = t + and is finite in the 
disk \z\ < 1, i.e. in the whole g 2 -plane. This ensures the correct analytic structure in the 
complex plane and the proper scaling, f+(q 2 ) ~ l/q 2 at large q 2 . 

As mentioned in the Introduction, f+(q 2 ) must satisfy also the condition lmf + (q 2 ) ~ 
(q 2 — £+) 3 / 2 near t + . Also, analyticity implies that near threshold Re f+(q 2 ) ~ a+ + b + (q 2 — 
t + ) + . . ., where a+ and b + are constants. We recall that from the definition (jlj) of the 
variable z = z(q 2 ,to) it follows that the threshold t + is mapped onto the point z = —1, and 
[z + 1) ~ const, x (q 2 — t + Y^ 2 near z = —1. Then, it is easy to see that /+ must satisfy the 
condition 



d^ 



= 0, (12) 
which, written in terms of the coefficients bk appearing in pip , takes the simple form 

K 

J2(-±) k+1 kb k (t )=0. (13) 
k=i 



By inserting in Eq. (fTTj) the solution of (fl3l) . written as 



we arrive at the expression 



,2 



K 

k=0 



Z k_(_ l) k-K» Z K 

K 



(15) 



where 2; = z(q ,t ). This is the parametrization that we investigate in the present work. 
As concerns the conformal mapping, i.e. the parameter t in ([111 , it was remarked in 
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321 ] that for t = £ op t with 

t opt = (m B + m 7r )(V^ - v 7 ^) 2 = 20.062 GeV 2 , (16) 

the semileptonic domain is mapped onto the symmetric interval \z\ < 0.279 in the z-plane. 
As we shall discuss below, this choice minimizes the maximum truncation error in the 
semileptonic domain. An additional argument for this choice is provided by the general 
theory of the representation of data distributed along an interval: as shown in [37], in this 
case the optimal expansion of the function is obtained by using a complete set of orthogonal 
polynomials of the variable that maps the original complex cut plane onto an ellipse, such 
that the cut becomes the boundary and the physical range is mapped onto the interval 
situated between the focal points. We checked that the optimal ellipse given in jsjjj is in our 
case very close to the circle \z(q 2 , t opt )\ = 1. 

Other choices of to are useful if one is interested in having a more accurate description 
in a specific energy range. Two values, to — and to = t- were investigated in [22j ]. In our 
analysis we shall adopt the choice t = t opt , with t opt given in ( TTTjl) . 



IV. UNITARITY CONSTRAINT 

The unitarity condition (J7j) can also be expressed in terms of the coefficients bk(to). By 
comparing the representations in Eqs. ([3]) and f[TT[) we have 

00 K 

J2dn(t )z n = ^(z)J2h(t )z k , (17) 

n=0 k=0 
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where ^f(z) is a known function 

9(z) = m% * *(z) (1 " Z)2{1 " Z * ? (18) 

We denote by z* = z(m 2 B ,,t ) the position of the pole in the variable z, and $(2) = 
0(g 2 (z),t o ) is the outer function expressed in terms of z, by using the inverse of Eq. (j4]) 

) 2 

1 - 

The function ^(2), which depends also on the parameter to, is analytic in \z\ < 1. Thus, we 
can expand it around z = as 

00 

= (20) 

fc=0 

Inserting this expansion in (ITTj) . we obtain 



g 2 (z)=t + -(t + -t ) (^) • (19) 



mia[ii",n] 

a n (*o) = X] Vn-k(to)bk(to), n > 0. (21) 

fc=0 

Then the inequality (jTJ), expressed in terms of the coefficients bj(t ), reads as 

K 

B jk (t )b,(t )b k (t ) < 1, (22) 

j,k=0 

where 

00 

B jk (to) = ^2r] n (t )r]n+\j~k\(to)- (23) 

n=0 

From Eqs. ([6]) and (fTSl) it follows that the function ^(z) is bounded in the closed disk \z\ < 1, 
so its Taylor coefficients rjj are rapidly decreasing. Therefore, the coefficients Bjk(t ) can be 
computed by performing in (1231) the summation upon n up to a finite order, about 100 in 
practice. 

As discussed in 



33| . the leading contributions to the sum over the coefficients a 2 , which 
appears in the unitarity condition ((7|), are of order (A/ 'nib) 2 i n the heavy 6-quark expansion. 
Thus, we expect that the 1 appearing on the right-hand sides of the constraints of Eqs. ([7]) 
and (1221) is a significant overestimate, the real bound being more realistically on the order 
of a few per mil. Were we to consider these more stringent bounds in the sequel, we would 
be able to reduce the number of terms kept in the series expansion of the form factor. 
This would significantly reduce the systematic uncertainty that we encounter when using 
this expansion to extrapolate the form factor to regions where it is not constrained by the 
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t (GeV 2 ) B 00 Bqi B 02 B 04: B 04: B 05 



0.0197 -0.0049 -0.0108 0.0057 0.0006 -0.0005 

t opt 0.0197 0.0042 -0.0109 -0.0059 -0.0002 0.0012 

t_ 0.0197 0.0118 -0.0015 -0.0078 -0.0077 -0.0053 



TABLE I: The matrix elements, Bjj-(to), which enter the unitarity bound (|22p for K = 5 and 
several values of to- The remaining coefficients are obtained from the relations £?j(j+fc) = Bq^ and 
the symmetry property Bj^ = B^j, obvious from Eq. (|23p . 

input that we use. However, in the absence of a more precise quantitative argument for 
strengthening the bound, we choose to keep Eq. (|22j) as it stands. In any event, if and when 
such an argument is found, the procedures explained in the following sections can be carried 
over as is, simply replacing the right-hand side of the inequality (|22|) by the relevant smaller 
number. 

For the numerical evaluation of the coefficients Bjk we need the value of Xi- (0) entering 

30] 



the outer function ([6]). Perturbative QCD and the operator product expansion give 
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_ 3[1 + l.Ua s (m b )} _ m b (uu) _ (a s G 2 ) 
Xl ~ [ } ~ 32vr 2 m 2 m 6 b Ylvm\ ' 1 ' 

where m& = 4.9 GeV is the pole mass and fh b the MS mass, with m&(2 GeV) ~ 4.98 GeV, 



obtained from mh(rh b ) ~ 4.2 GeV |38j and the four-loop running in the MS scheme [39]. 
We took a s (rhb) = 0.22 (3^] - The gluon condensate has the standard value (a s G 2 ) = 
0.038 GeV 4 given in [40|], while for the quark condensate we used the two-flavor value (uu) ~ 



— (278 MeV) 3 in the MS scheme at scale 2 GeV 41(. From the above values we derive 
rhb (uu) ~ —0.107 GeV 4 at scale 2 GeV, and we adopt this value also at scale m^, since the 
scale dependent corrections to the product are negligible. By inserting in ( |24l) the above 
central values we obtain ~ 5-01 x 10~ 4 . For illustration, we give in Table U the 

coefficients Bjk(to) calculated with this value of Xi-(0) f° r K = 5 and several values of to- 

We mention also that when K —>■ 00, the expansion in Eq. ffTTJ) is convergent in the whole 
disk \z\ < 1, i.e. in the whole g 2 -plane cut along the real axis for q 2 > t+. Moreover, the 
unitarity condition fl22l can be used to derive an explicit upper bound on the truncation 
error. We present the derivation of this bound in the Appendix. 
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V. THEORETICAL AND EXPERIMENTAL INPUT 



At low q 2 , the form factor is calculated in the frame of LCSR 



3, 



11 



121. We use 



/lcsr = /+(0) = 0.26, with the uncertainty <5/lcsr — 0.03 [22J. Lattice calculations provide 
;he value of the form factor at eight additional g 2 -points: three are taken from FNAL-MILC 
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, 3jJ and five from HPQCD, updated in [17]. As in 2, 4], we take the three FNAL-MILC 
results from [l]. 1 

The available experimental data consist in the partial branching fractions over bins in q 2 . 
We use 10 data from the tagged analyses (4 bins from CLEO [6|, which replace the older 
data |||, 3 from Belle and 3 from BaBar [8]), and 12 bins from the untagged BaBar 
analysis [9], where the full covariance matrix is available. The total number of data points 
from theory and experiment is 31. 

It is convenient to define the global % 2 



X 2 (h,\Vub\) = Xth + X 



exp 



(25) 



where 



Xth 



X, 



exp 



£ [fT - Uif$)\C;i\ff - /+(?*)] + (/+(«) - /lcsr) 2 /(^lcsr) 2 , 
j,k=l 

22 

= E - B AU)\C-B jk \Bt - B k {U)\ . (26) 
j,k=i 



In the above relations, the /- n denote the values of the form factor calculated on the lattice 
at the points q 2 ; BJ 1 are the experimental partial branching fractions and Bj(f + ) the values 
calculated by integrating Eq. (j2J) over the bins [q 2 ,q 2 +1 ], with a given parametrization for 
the form factor f + (q 2 ). To convert to a rate, we use the B° lifetime, = l/T tot = (1.527 ± 
0.008) x 10~ 12 s [42]. 



The covariance matrices C and Cb are written formally: in practice they are block 
diagonal, with independent blocks for each independent set of experimental data or lattice 



1 As we were finalizing this work, Fermilab and MILC presented, in [18], a substantial update of their 
lattice calculation of the form factor f+{q 2 )- Since their new results agree within errors with those of 
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161 ] and since our goal here is to illustrate the workings of our new parametrization, we have chosen 



not to update our analysis. Instead, we encourage the authors of [18| to perform their analysis with our 
improved parametrization. 
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results. Unfortunately, the covariance matrices are not provided for the lattice calculations. 
Thus, we make here a set of reasonable assumptions on the possible correlations based on 
the information provided in t he p apers and on our experience with such calculations. The 



lattice results of 



16| and 17J are obtained using different discretizations for the heavy 
quark, but on subsets of the MILC, Nf = 2 + 1, gauge configurations which have significant 
overlap. Thus, in addition to assuming that the statistical errors on f+{q 2 ) at different q 2 
within each calculation have a 50% correlation, 2 we assume that there is a 25% correlation 
between the errors in the two calculations. Such correlations in the statistical errors have 
been assumed to be negligible in previous work . Regarding the systematic errors, 

since the heavy-quark discretizations and methods used are different, we assume negligible 
correlations in the systematic errors between the two calculations, but given their nature, 
assume 100% correlations within each simulation. Though these assumptions cannot replace 
covariance matrices determined by the lattice collaborations themselves, we believe that they 
are reasonable and will not lead to underestimated errors. We have verified that without 
these correlations, for instance, the results for / + versus q 2 quoted below would have fit 
errors reduced by up to 30%. 



VI. RESULTS OF THE FITS 

We performed a combined fit of the above input by minimizing \ 2 defined in ( 1251) . using 
the representation of f+(q 2 ) given in (1T51) . with z = z(q 2 ,t opt ). The free parameters are \ V u b\ 
and the real coefficients bk, k < K — 1 subject to the unitarity constraint (1221) . with bx given 
by the expression ()14p . The total number of par ameters is iV = K + 1. 



According to convex optimization theory 



43| . the optimum values b% and the optimal 



Lagrange multiplier Ao minimizing the Lagrangian 
C(b J ,\V ub \) = X 2 



(bj, \V ub \) + B ^ bk ~ 1 ) ' ( 2? ) 



2 Points at different q 2 within a given simulation are obtained on the same statistical ensemble with very 
similar methods and are thus expected to be strongly correlated. A glance at Fig. [2 in which the lattice 
results are plotted, should convince the reader that such correlations are present. 
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q [GeVI 



10 15 
q 2 [GeV 2 ] 



FIG. 1: Left: the form factor f+{q 2 ) computed with the representation (|15p and the parameters 
from (I29p -(|3ip. Right: the form factor calculated with the traditional parametrization ([3]) and the 
parameters from (|32p-(j34p. 



satisfy the alignment condition 



A (fB,4°'6f-l)=0. 

\j,k=0 J 



(28) 



Therefore, either Ao = and the optimal parameters b^ of the unconstrained minimum of 
X 2 (bk, \Vub\) satisfy automatically the constraint ( 1221) . or Ao 7^ 0, when the optimal parame- 
ters saturate the constraint fl22l) . 

A nontrivial form factor is obtained for K > 2, i.e. a total number of parameters, N > 3. 
The results of the fits obtained by increasing N are presented below: 

N = 3; K = 2; X 2 = 16.87, x7 dof = 0.60 
b = 0.437 ± 0.024, bt = -0.552 ± 0.066, 
\V ub \ = (3.47 ±0.21) x 10~ 3 . 



(29) 



N = 4; K = 3; x 2 = 16.80, x7 dof = °-62 
6 = 0.438 ± 0.024, 6i = -0.496 ± 0.234, b 2 
\V ub \ = (3.45 ±0.23) x 10~ 3 



-0.459 ± 0.733 



(30) 



N = 5; K = 4; X 2 = 16.11, X 2 / dof = 0.62 

b = 0.444 ±0.025, 61 = -0.439 ± 0.249, b 2 = -2.111 63 = 3.928^;^ 

\V ub \ = (3.50 ± 0.24) x 10~ 3 . (31) 
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For simplicity we omitted the upper index "(0)" in the notation of the optimal parameters. 
All the errors indicated are statistical. We mention that for the fits (1291) and (130]) the 
unitarity constraint (I22p is not saturated, while the values in (131 p saturate the inequality 
(|22p . The asymmetric errors on the coefficients b 2 and 63 in f[3"Tl) are produced mainly by 
this constraint. 

The form factor calculated with the central values of the parameters from (!29l) - (l3~T|) 
is plotted in the left panel of Fig. [TJ For comparison we repeat the analysis also with 
the standard parametrization ([3]). We recall that this parametrization has an unphysical 
singularity at threshold, and we cannot impose the threshold condition that we use above. 
Therefore, for a certain K there are K + 1 parameters a^, constrained by the unitarity 
condition ([7]), and the total number of parameters is N = K + 2. The best fits for the lowest 
values of N are 

N = 3; K = 1; \ 2 = 19.33, xV^oi = 0.69 
a = 0.025 ± 0.0041, a x = -0.039 ± 0.0042 

\V ub \ = (3.53 ± 0.22) x 10^ 3 (32) 

N = 4; K = 2; X 2 = 16.97, x7 d o f = 0.63 

a = 0.026 ± .001, 01 = -0.021 ± 0.014, a 2 = -0.076 ± 0.053 

\V ub \ = (3.43 ± 0.23) x 10~ 3 (33) 
N = 5; K = 3; x 2 = 16.01, x7 d o f = 0.62 

a = 0.027 ± 0.002, a x = -0.019 ± 0.014, a 2 = -0.195 ± 0.139, a 3 = 0.409 ± 0.421 
= (3.51 ± 0.25) x 10' 3 . (34) 

In all cases, the unitarity constraint (171) is not saturated. The corresponding form factor 
is plotted in the right panel of Fig{TJ For large values of q 2 the results indicate a more 
pronounced variation with iV than that of the curves in the left-hand panel of the figure. 



VII. SYSTEMATIC ERROR 



In the present framework, the systematic error on the values of the form factor must 
account for the effect of truncating the expansion (fTTl) at a finite order K. As shown in the 
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appendix, the unitarity constraint (|22|) can be exploited to derive an upper bound on the 
truncation error. However, this estimate is too conservative for low values of K. A more 
realistic prescription is given by the magnitude of the next term in the expansion, allowed by 
the unitarity constraint. Denote by the maximum value of the modulus allowed 
by the condition (I22p . for fixed values of k < K, given by the fit. We note that, although 
the inequality (|22|) may be saturated by the latter values, as happens with the values in Eq. 
fl3Tj) . a nonzero value for is obtained, since the convex condition (122]) is not a sum of 
squares. 

According to the above discussion, we adopt as a realistic prescription for the systematic 
error on the form factor the quantity 

Sf+(q 2 )syst = 7^7-2 > ( 35 ) 

where z = z(q 2 ,t opt ). Using the optimal b^ from (|2"91) - (|3~TI) . bx from (TBI) and the values of 
Bjk for t = topt given in Table IH we obtain 

b max = 7 b max = j ^ ^ax = 3 Q2 ^g^ 

With these values, the numerator of fl3"oT) calculated at the limits of the physical region, 
\z m ax\ = 0.279, is a fraction of 35.9%, 10.3% and 1.2% from the corresponding first coefficient 
bo given in (!29"]) - (!3~T|) . for A^ = 3, A^ = 4 and A^ = 5, respectively. For illustration we give 
also the values of the form factor and errors at the highest point t_, where the systematic 
error defined in (1351) has the largest value 

f + (t_) = 8.34 ± 0.53 (stat) ± 2.30 (syst), N = 3, 
f + (t_) = 7.96 ± 1.61 (stat) ± 0.66 (syst), N = A, 
f + (t_) = 4.97 tlf 3 (stat) ± 0.07 (syst), N = 5. (37) 

For N = 5 the systematic error is very small. Actually, it is negligibly small compared to 
the statistical error along the whole physical region. By going up to A^ = 5, we can neglect 
the systematic error altogether for the determination of V u b and for the form factor in the 
physical region. We shall adopt this choice as our optimal parametrization. 
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VIII. BEST PARAMETRIZATION IN THE PHYSICAL RANGE 



As discussed above, we adopt the expansion (fT5|) for K = 4, which writes as 

1 3 

/+o? 2 ) = -, — 2/ 2 y2b k 

l-q 2 /m z B , 

The best parameters and their statistical errors, already given fl3T|) . are 

6 = 0.444 ± 0.025, h = -0.439 ± 0.249, b 2 = -2.111 ±g;$g, 6 3 = 3.928 (39) 

We mention that the parameters are correlated, the correlations being highly nongaussian 
because of the unitarity constraint. 

As shown in ( I3T1) . the fit gives x 2 = 16.11 and x /dof = 0.62. For completeness we list 
below the separate contributions to Xth an d xtxp °f the various data sets, compared with 
the number n of points: 

Xlcsr = 0.032 (n = 1) Xfnal-milc&hpqcd 
xLiie = 0-017 (n = 3) Xcleo = 2.42 (n = 

XBaBar-t = 4 -04 (/l = 3) xlaBar-u = 8 - 93 

The description of all the sets is very good, except for the BaBar tagged (t) data |8j, where 
X 2 is larger than the number of points. 

The form factor calculated using the expression (138]) and the parameters from Eq. fl39|) 
is shown in Fig. [21 where in the right panel we plot the z polynomial in the numerator of 
Eq. (138]) . The error bands represent the purely statistical error. We emphasize that we 
do not use the linear approximation in the error propagation, but apply the standard A% 2 
analysis, i.e. by finding the range of variation of a given parameter corresponding to a 
change in x 2 , minimized over all other parameters, by one unit. The unitarity constraint 
plays a nontrivial role, being responsible for the asymmetric errors, especially near the right 
end of the semileptonic range. For completeness, the values of the form factor are given in 
Table HT1 for a sample of q 2 in the semileptonic domain. In Table HTT1 we compare the results 
of our combined fit with the lattice results used as input. 

As seen from the values given in (1371) . the gradual reduction of the systematic error with 
the increase of iV is balanced by the increase of the statistical error. It is of interest to 
compare the total error on the values of the form factor for N = 5 and N = 4. In Fig. [3] 
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z k - (-l) k -z 4 



(38) 



= 0.68 (n = 3 + 5) 
= 3) 

= 12). (40) 
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FIG. 2: Left: the form factor / + (g 2 ) given by Eqs. (|38p and (|39p . The error band corresponds 
to statistical: in our approach, systematic errors coming from the parametrization are negligible 
in the semileptonic domain, as explained in the text. The theoretical LCSR result from [12] and 
the lattice results from [la . \l\ are also shown. Right: the numerator in (|38p for the optimal 
parametrization. 
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FIG. 3: Variation of the error in the physical region: solid/dashed line indicates the difference 
between the plus/minus error for N = 5 and the error for N = 4. 

we plot the difference between these two errors as a function of q 2 . For N = 5 the error 
is purely statistical, for N = 4 it is calculated by adding quadratically the statistical and 
the systematic errors, the later one obtained from (1351) . As we already noted, the error 
in the N = 5 case is not symmetric due to the unitarity constraint, therefore we present 
separately the difference between the "plus" and "minus" N = 5 errors and the N = 4 
one. Figure [3] shows that the difference between the errors is practically zero for most of 
the energy range, including the energies where lattice input is available. At low values of 
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g 2 (GeV 2 ) f + {q 2 ) g 2 (GeV 2 ) f + (q 2 ) 



0. 


0.265 


+0.019 
-0.019 


18. 


1.142 


+0.071 
-0.070 


2. 


0.292 


+0.019 
-0.019 


19. 


1.304 


+0.075 
-0.076 


4. 


0.327 


+0.024 
-0.022 


20. 


1.503 


+0.083 
-0.084 


6. 


0.370 


+0.029 
-0.025 


21. 


1.749 


+0.098 
-0.076 


8. 


0.426 


+0.036 
-0.031 


22. 


2.059 


+0.132 
-0.131 


ID 


498 


+0.043 
-0.036 




2 453 


+0.223 
-0.200 


12. 


0.593 


+.0490 
-0.044 


24. 


2.967 


+0.456 
-0.337 


14. 


0.720 


+0.056 
-0.052 


25. 


3.644 


+1.033 
-0.611 


16. 


0.895 


+0.063 
-0.060 


26. 


4.535 


+2.579 
-1.226 


17. 


1.007 


+0.067 
-0.065 


26.42 


4.967 


+3.973 
-1.726 



TABLE II: The form factor at a variety of q 2 values in the semileptonic domain, as obtained 
with the expression of Eq. ([38]) and the parameters given in Eq. ([39]) . The errors are obtained as 
described in the text. 

q 2 , in particular at q 2 = 0, the total error decreases when passing from N = 4 to N = 5. 
On the other hand, at high values of q 2 the plus error for iV = 5 is larger than the N = 4 
error. Qualitatively, the fact that the unitarity constraint increases the error band in the 
upward direction for iV = 5 is in agreement with the shape of the form factor for N = 3 and 
N = 4, seen in the left panel of Fig. [TJ One may ask whether it is not preferable to take as 
a best prediction the parametrization with N = 4. In our opinion, this is not the case: the 
advantage of our prescription is that the systematic errors are negligible along the whole 
physical region. Thus, we avoid any bias related to the specific form of the truncation error 
for the determination of the form factor and of V u b. Our results show that a representation 
of the form factor having a small uncertainty over the whole physical region, including its 
upper end, is not possible with the present input information. 

From the above comment we expect even larger errors if the expression ( [38]) is used to 
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g 2 (GeV 2 ) / + (g 2 )i a ttic e /+(g 2 )fit 



15.87 0.799 ± 0.058 ± 0.088 0.882 ±0.059 



18.58 1.128 ± 0.086 ± 0.124 1.232 



+0.073 
-0.074 



24.09 3.263 ± 0.324 ± 0.359 3.020 ±0.354 



17.34 1.101 ± 0.053 ± 0.105 1.050 tg-jgg 



18.39 1.273 ± 0.099 ± 0.121 1.201 



+0.072 
-0.072 



19.45 1.458 ± 0.142 ± 0.139 1.388 tg^re 



20.51 1.627 ± 0.185 ± 0.155 1.622 ±o.olo 



21.56 1.816 ±0.126 ± 0.173 1.914 ±o*.ii3 



TABLE III: Comparison of the form factor / + (g 2 ) obtained at a variety of recoils in the lattice 
computations of FNAL-MILC Q, Q (rows 1-3) and HPQCD Q (rows 4-8), with the results of 
the combined fit in Eq. (|39p . to the parametrization given in Eq. (|38p . In the lattice results, the 
first error is the combined statistical and chiral extrapolation error, while the second is an 11% 
systematic for FNAL-MILC Q and a 9.5% error for HPQCD[17[. 



calculate the form factor outside the physical region. In particular, we consider the residue 
of f+(q 2 ) at the pole q 2 = m 2 Bt) denned as 

r + = lim (l-q 2 /m 2 B ,)U(q 2 ). (41) 

q 2 — >m,g t 

The systematic error, calculated using the prescription fl35|) . is no longer negligible at the 
position = —0.504 of the pole. From ([38"]) and (|39|) we obtain 

r + = -0.12 tof 4 (stat) ± 0.10 (syst). (42) 

The large statistical and systematic errors show that a reliable extraction of the residue from 
the extrapolation of our best fit is not possible. Additional information on the behavior of 
f+(q 2 ) outside the physical region, like the absence of zeros, expected on general grounds for 



form factors [44j], or the monotony, valid in some models 45), might improve the prediction. 
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Alternatively, one can use a different conformal mapping, i.e. a different value of the 
parameter to in which allows a better accuracy in the high energy range. A reasonable 
choice is to — t-, when the physical region (0, t_) is mapped onto the interval (0.518, 0) 
of the z plane, and the pole becomes z(m 2 B *,t_) = —0.262. Since the pole is closer to the 
origin, the systematic error at this point is now negligible already for the best fit with N = 4 
parameters, when we obtain 

r+ = 0.52 ±0.34 (stat). (43) 
This value can be converted to a prediction for gB*Bn = 2r + rriB* / /#* ■ Using, for instance, 



fa* = 0.196 ± 0.031 GeV [46j, we obtain 



g B . B n = 28.23 ± 18.46 (stat) ± 5.7 (6f B *), (44) 



to be compared with the lattice result Qb*Btt = 47 ± 3 (stat) ± 9 (syst) from [471 ] . A more 
detailed study of this problem will be presented in a future work. 

Before ending this section, let us make a few more comments on the standard analytic 
parametrization ([3]). We presented the results of fits to this parametrization in (l32l) - (!34|) 
and in Fig. [IJ As discussed in Sec. [Til the parametrization ([3]) has a fake singularity at the 
unitarity threshold, q 2 = t + , which is expected to produce distortions in the behavior of the 
form factor at large values of q 2 . For instance, from the fit with N = 5 parameters, i.e. 
using four terms in the expansion ([3]) and the best values from (134"1) . we obtain 



= 2.47 (stat) ± 1.69 (syst), [Eqs.©, (EI)], (45) 

and the residue 

r + = -2.18^-88 ( s t a t) ± 1.91 (syst), [Eqs.©, (jM])]. (46) 

The larger systematic errors are explained in part by the fact that now the expansion has 
only four terms (unlike in ( l38i) . where an additional term was introduced using the threshold 
condition). The statistical errors are also larger than for the iV = 5 fit which uses our new 
parametrization, Eqs. (137|) and fj4"2]) . showing that the singularity at threshold affects the 
behavior of the form factor near this point. 
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IX. DETERMINATION OF \V ub \ 



As shown in Sec. IVI} \V u b\ is one of the parameters determined by our fit: the optimal 
value and the statistical error are given in Eq. f[3"Tj) . We chose the parametrization such 
as the systematic error can be neglected along the whole physical region. Therefore, the 
determination of \ V u b\ will be free of systematic uncertainties. Adding an experimental error 
of 0.01 x 10~ 3 associated with the uncertainty in the B° lifetime [421 ]. our final prediction is 

\V ub \ = (3.50 ± 0.24) x 10" 3 . (47) 

This result depends of course on the theoretical and experimental input used, and will 
become more and more accurate as this input will improve. Our purpose in this work was 
mainly to prove the advantages of the simple parametrization Eq. (1381) of the form factor, 
which we recommend as a useful tool in future data analyses. 

The result (1471) is consistent with the most recent prediction from exclusive B — > ir decays 
4j. However, as discussed in Sec. Ull the analysis in [4| is based on a parametrization that 
does not fully satisfy the constraints of analyticity and unitarity. Moreover, the statistical 
correlations in the lattice results are neglected there. Thus, our analysis puts the extraction 
of \V u b\ from exclusive B —>■ n decays on a more rigorous basis. 

From the fit given in ( 13~TT) we obtain also 

|Kft|/+(0) = (9.31 ± 0.60) x 10" 4 , (48) 
to be compared with the result |V^b|/ + (0) = (7.6 ± 1.9) x 10~ 4 obtained with SCET and 



factorization 



48|. 



X. CONCLUSIONS 



We proposed a simple analytic parametrization for the semileptonic B — > tx vector form 
factor f + (q 2 ), by multiplying the factor accounting for the B* pole with a convergent ex- 
pansion in powers of a conformal mapping variable. The parametrization has the correct 
behavior at the unitarity threshold and satisfies perturbative scaling and the constraint de- 
rived from the positivity of the correlation function of the u'j^b current and its Hermitian 
conjugate. The latter was used also to define the systematic error due to the truncation 
of the expansion. By increasing up to K = 4 the number of terms in the expansion, we 
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obtained the representation given in Eqs. (!38|) - (l39|) . where the systematic error can be ne- 
glected along the whole physical region. From the combined fit of our parametrization to 
experimental results for the differential decay rate and to theoretical results for the form fac- 
tor, we obtained a prediction for \V u i\ given in (|47p . Our result confirms that \V u b\ extracted 
from the exclusive B — > tx decays is consistent with the global fits of the CKM matrix {49) ]. 
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APPENDIX 

In this Appendix we show that unitarity allows one to derive a bound on the remainder 
of the expansion (Till , defined as 



For simplicity we omit the dependence of the coefficients and the variable z upon t 0} which 
is kept fixed in the expressions given below. 

Using (j3J) and (iTTj) we express each coefficient bk as 




(A.l) 



k 



j=0 



k > 0, 



(A.2) 



where fjj appear in the expansion 



00 




(A.3) 



fc=0 



with ^(z) given in Eq. (Tl8i) . By the Cauchy inequality we obtain from (1A.2I) 



k 




(A.4) 



j=0 



3=0 
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and using (J7J) 

k 

|& fc | < {]>>f} 1/2 , k >°- ( A - 5 ) 

3=0 

Therefore, the remainder (1A.1I) is bounded in terms of calculable quantities 

1 oo k 

The upper bound (IA.6[) can be made sufficiently small for a certain K and |^| < 1. This 
follows from the properties of the coefficients fjj defined in flA.31) : indeed, the function 
l/^>(z) has singularities on the boundary \z\ = 1, but it is analytic inside the disk \z\ < 1. 
Therefore, although the Taylor coefficients fjj increase with j, the increase is such that sum 
Ylk>KVk\z\ h can be made arbitrarily small for a certain K and \z\ < 1. The same is true 
for the coefficients appearing in ( 1A.6I) . this fact being obvious, in particular, if we use the 

k 

upper bound {Y^ Vj} 1 ^ 2 < k Vk, valid for sufficiently large k. Using this estimate, the sum 

3=0 

in flA.61) is related to the remainder of the Taylor expansion of the derivative of the function 
l/ty(z), which can be made arbitrarily small since the series is convergent for \z\ < 1. 
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